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Internal modes of discrete solitons near the 
anti-continuum limit of the dNLS equation 

Dmitry Pelinovsky and Anton Sakovich 
Department of Mathematics, McMaster University, Hamilton ON, Canada, 

Abstract 



Discrete solitons of the discrete nonlinear Schrodinger (dNLS) equation are compactly sup- 
ported in the ant i- continuum limit of the zero coupling between lattice sites. Eigenvalues of the 
l/-v ' linearization of the dNLS equation at the discrete soliton determine its spectral stability. Small 

eigenvalues bifurcating from the zero eigenvalue near the anti-continuum limit were character- 
ized earlier for this model. Here we analyze the resolvent operator and prove that it is bounded 
in the neighborhood of the continuous spectrum if the discrete soliton is simply connected in 
the anti-continuum limit. This result rules out existence of internal modes (neutrally stable 

c^ . eigenvalues of the discrete spectrum) near the anti-continuum limit. 
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1 Introduction 

The discrete nonlinear Schrodinger (dNLS) equation is a mathematical model of many physical 
J^ phenomena including the Bose-Einstein condensation in optical lattices, propagation of optical 

pulses in coupled waveguide arrays, and oscillations of molecules in DNAs |12j . Discrete solitons 
(stationary localized solutions) are used to interpret the results of physical experiments and to 
characterize global dynamics of the dNLS equation with decaying initial data. 

Discrete solitons are compactly supported in the anti- continuum limit of the zero coupling be- 
tween lattice sites. Different families of discrete solitons can be uniquely characterized near the 
anti-continuum limit from a number of limiting configurations pQ . This is the main reason why the 
anti-continuum limit has been studied in many details after the pioneer works of Eilbeck at al. |10] 
and Aubry & Abramovici [2]. The existence of discrete solitons (also called discrete breathers in 
the context of the discrete Klein-Gordon equation) was rigorously justified with implicit function 
theorem arguments by MacKay & Aubry |15j . Their work on existence of discrete solitons led to fur- 
ther progress in understanding their stability properties as well as nonlinear dynamics of nonlinear 
lattices [3|0. 

Spectral stability of discrete solitons is determined by eigenvalues of the discrete spectrum of 
an associated linearized operator because its continuous spectrum is neutrally stable. Unstable 
eigenvalues can be fully characterized near the anti-continuum limit because they bifurcate from the 
zero eigenvalue of finite multiplicity and the zero eigenvalue is isolated from the continuous spectrum. 
Characterization of unstable eigenvalues for each family of discrete solitons bifurcating from a 
compact limiting solution was obtained by Pelinovsky et al. |17j with an application of Lyapunov- 
Schmidt reduction technique. Beside the unstable eigenvalues, the same technique was used to 
characterize a number of neutrally stable eigenvalues of negative energy (also called eigenvalues of 
negative Krein signature) which bifurcate from the same zero eigenvalue. These isolated eigenvalues 
of negative energy may become unstable far from the anti-continuum limit because of collisions with 
eigenvalues of positive energy (also called internal modes) or with the continuous spectrum of the 



linearized operator. Isolated eigenvalues of negative energy may also induce nonlinear instability if 
their multiples belong to the continuous spectrum [7]. 

In the same anti-continuum limit, another bifurcation occurs beyond the applicability of the 
Lyapunov-Schmidt reduction technique: a pair of semi-simple nonzero eigenvalues of infinite mul- 
tiplicity transforms into a pair of continuous spectral bands of small width. This transformation 
may produce additional eigenvalues of the discrete spectrum similar to what happens for the dis- 
crete kinks (which are non-compact solutions of the nonlinear lattice in the anti-continuum limit) 
|18| . No complex unstable eigenvalues may bifurcate from the semi-simple nonzero eigenvalues of 
infinite multiplicity because such eigenvalues are excluded by the count of unstable eigenvalues in 
|17j . Nevertheless, internal modes may in general be expected outside the continuous spectrum. 

It is important to know the details on existence of internal modes because of several reasons. 
First, these internal modes may collide with eigenvalues of negative energy to produce the Hamilton- 
Hopf instability bifurcations [17J. Second, analysis of asymptotic stability of discrete solitons de- 
pends on the number and location of the internal modes [9j [13]. Third, the presence of internal 
modes may result in long-term quasi-periodic oscillations of discrete solitons [8|. 

In this paper, we address bifurcations of internal modes from semi-simple nonzero eigenvalues of 
infinite multiplicity. We continue the resolvent operator across the continuous spectrum and prove 
that it is bounded near the end points of the continuous spectrum if the discrete soliton is simply 
connected in the anti-continuum limit, see Definition [2j As a result, no internal modes exist in the 
neighborhood of the continuous spectrum. These results hold for any discrete soliton of the dNLS 
equation with any power nonlinearity near the anti-continuum limit. 

There are multiple numerical evidences that no internal modes exist near the anti-continuum 
limit for the fundamental discrete soliton, which is supported at a single lattice site in the zero 
coupling limit. In particular, this fact is suggested by Figure 1 in Johansson &: Aubry |11| and by 
Figure 2.5 in Kevrekidis |12j . Our article presents the first analytical proof of this phenomenon. 

The paper is organized as follows. Section 2 reviews results on existence and stability of discrete 
solitons near the anti-continuum limit. Section 3 is devoted to analysis of the resolvent operator with 
the limiting compact potentials. Section 4 develops perturbative arguments for the full resolvent 
operator. Section 5 considers a case study for the resolvent operator associated with a non-simply- 
connected 2-site discrete soliton. Appendix A is devoted to the cubic dNLS equation, for which 
perturbation arguments are more delicate. 

Notations. We denote the bi-infinite sequence {u n } n £z by u. The l p space for sequences is 
denoted by P(Z) and is equipped with the norm 
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The algebraically weighted space Zf (Z) with sGlis the P(Z) space for the sequence 

{(l+n 2 y/ 2 u n } nez . 
A disk of radius 5 > centered at the point Ao S C on the complex plane is denoted by B$(\q) C C. 

2 Review of results on discrete solitons 

Consider the dNLS equation in the form 

iii n + e(n n+ i - 2u n + u„_i) + \u n \ 2p u n = 0, n G Z, (1) 



where the dot denotes differentiation in t G R, {u n (t)} n <=z '■ R Z — ?• C is the set of amplitude functions, 
and parameters e G R and p G N define the coupling constant and the power of nonlinearity. The 
anti-continuum limit corresponds to e = 0, in which case the dNLS equation JT]) becomes an infinite 
system of uncoupled differential equations. 

Discrete solitons are defined in the form u n (t) = 4> n e lt , where the frequency is normalized thanks 
to the scaling symmetry of the power nonlinearity. By the standard arguments [16] based on the 
conserved quantity 

e / : 4>n4>n+i - 4>n4>n+i = const in neZ, (2) 

it is known that if {4> n }neZ decays to zero as \n\ — > oo, then {4> n }neZ is real-valued module to 
multiplication by e td for any 9 G R. The real-valued stationary solutions are found from the second- 
order difference equation 

(1 - <$> n = e(<£ n+1 - 2</> n + </>„_!), n G Z. (3) 

The algebraic system is uncoupled if e = 0. 

Let us consider solutions of the difference equation ([3]) for cf> G / 2 (Z). If e = and pgN, the 
limiting configuration of the discrete soliton is given by the compact solution 



e = : 0W = ^ Sn - ^ S n , (4) 

ng£/+ neC/_ 



where U± are compact subset of Z such that {/+ D U- = and S n is the standard unit vector in 
/ 2 (Z) expressed via the Kronecker symbol by 

\Pn)m = On,rm m G IL. 

We will denote the number of sites in U± by \U±\. The following proposition gives a unique analytic 
continuation of the compact limiting solution (|5J) to a particular family of discrete solitons (see 
[El CE1I7] for the proof). 



Proposition 1 Fix U + , U- C Z such that U + P\U- = and \U+\ + \U-\ < oo. There exists eo > 
such that the stationary dNLS equation ([3J) with e G (— eo,£o) admits a unique solution cf) G / 2 (Z) 
near (jy- ' G / 2 (Z). T/ie map (— eo,eo) 9 e i— >• 4> G Z 2 (Z) is analytic and 

3C>0: ||0-0(°)[|,a<C|e|. (5) 

Moreover, there are k > and C > suc/i i/toi /or any e G (— eo, eo) 

|<£n| <Ce- K|n| , nGZ. (6) 

Remark 1 Thanks to the exponential decay ^j, the solution <f> G Z 2 (Z) of Proposition^ belongs to 
4> G / 2 (Z) /or any s > 0. 

By Proposition [TJ the solution for a given 0' ' can be expanded in the power series 

oo 

= 0(°)+JV0« eG(-e ,e ), (7) 

fe=i 

where correction terms {</r '}keN are uniquely defined by a recursion formula. 



Spectral stability of the discrete solitons is determined from analysis of the spectral problem 

L_|_u = —Aw, L_w = Au, (8) 

where A G C is the spectral parameter, (u, w) G l 2 (Z) x Z 2 (Z) is an eigenvector, and L+ are discrete 
Schrodinger operators given by 



(L+u) n = -e(u n+1 - 2u n + u n _i) + (1 - (2p + l)4>l p )u n 
(L_w) n = -e(w n+ i - 2w n + w n -i) + (1 - <thF)w n , 



2 n G Z. (9) 



We recall basic definitions and results from the stability analysis of the spectral problem (jSJ). 

Definition 1 XTte eigenvalues of the spectral problem (0|) with Re(A) > (resp. Re(A) = 0) are 
called unstable (resp. neutrally stable). If A G iR. is a simple isolated eigenvalue, then the eigenvalue 
A is said to have a positive energy if (L+u, u)p > and a negative energy if (L+u, u)p < 0. 

Remark 2 If X £ iM is an isolated eigenvalue and (L+u, u)p = 0, then A is not a simple eigen- 
value. In this case, the concept of eigenvalues of positive and negative energies is defined by the 
diagonalization of the quadratic form (L+u, u}/2, where u belongs to the subspace ofl 2 (Z) associated 
to the eigenvalue A of the spectral problem |5|) and invariant under the action of the corresponding 
linearized operator (see \6\ for the relevant theory). 

In the anti-continuum limit e = 0, the spectrum of L+ (resp. L_) includes a semi-simple 
eigenvalue — 2p (resp. 0) of multiplicity N = |L+| + |E/_| < oo and a semi-simple eigenvalue 1 of 
multiplicity |Z\{[/+ U U-}\ = oo. The spectral problem ([8]) has a pair of eigenvalues A = ±i of 
infinite multiplicity and the eigenvalue A = of geometric multiplicity N and algebraic multiplicity 
2N. The following proposition describes the splitting of the zero eigenvalue near the anti-continuum 
limit for e > (see |17] for the proof). 

Proposition 2 Fix U+, [/_CZ such that U + n U- = and N := \U + \ + \U-\ < oo. Fix e > 
sufficiently small and denote the number of sign differences of {(j)n }n£U+uU- by uq. 

• There are exactly no negative and N — 1 — hq small positive eigenvalues of L_ counting 
multiplicities and a simple zero eigenvalue. 

• There are exactly uq pairs of small eigenvalues A £ iM. and N—1 — uq pairs of small eigenvalues 
A G R of the spectral problem |||) counting multiplicities and a double zero eigenvalue. 

Proposition [2] completes the characterization of unstable eigenvalues and neutrally stable eigen- 
values of negative energy from negative eigenvalues of L + and L_ . In particular, we know from [6j 
that if Ker(L + ) = {0}, Ker(L_) = span{</>}, and (L+ (f>,(f))i2 ^ 0, then 

n(L + )- Po = N- + Nr+N c , 

n(L_) = N+ + Nr + N c ' ' 



»C) 



where n(L±) denotes the number of negative eigenvalues of L±, Nf denotes the number of eigen- 
values A G iM with negative energy, N c denotes the number of eigenvalues with Re(A) > and 
Im(A) > 0, N^ 1 denotes the number of eigenvalues A G R with (L+u, u)p ^ 0, and 

r i if (l + v,0)p<o, 
po \ o if (l;V,0)p>o. 



To compute po, we extend the family of discrete solitons by parameter w as solutions of 
(w - (fil p )4> n = e((j) n+1 -2<t> n + (j) n -i), n£Z. 
Differentiation of equation (jlip in u) at u) = 1 gives 



(11) 



(l;V,0)p = -(9 ( 



•u<P\u=l 



,<t>)p 



l_d_ 

2dw 



II0IIP 



N 



+ 0(e), 



where in the last equality we used Proposition Q] and the anti-continuum limit 

e = 0: \\4>{u)\\l=Nu l ' p . 

Therefore, po = 1 for small e > 0. 

By Proposition [2j we have n(L_) = no and Nf > no- Also, n{L + ) = N . Using the count {[TO 
we have for small e > 



AC 



0, NT = N-1- n , NT = n , A c = 0. 



(12) 



Equality (|12p shows that besides the small and zero eigenvalues described by Proposition [2j 
the spectral problem may only have the continuous spectrum and the eigenvalues on iM. with 
positive energy. These eigenvalues of positive energy are called the internal modes and existence of 
such eigenvalues for small e > is the main theme of this article. 

3 The resolvent operator for the limiting configuration 

Let us consider the truncated spectral problem ([8]) after is replaced by (fy-K The resolvent 
operator is defined from the inhomogeneous system 



-e(u n+1 - 2u n + u n -\) + u n - (2p + 1) J2 



m£U+UU- 



On,mUm T AW n — r n) 



-e(w n+ i - 2w n + Wn-l) +W n - Y^mdU+YJU- dn,mW m ~ A"U n = G n , 



n e N, (13) 



where F, G G l 2 (^) are given. Since we are interested in the continuous spectrum and eigenvalues 
on iM., we set A = — iCl and use new coordinates 



a n := u n + iw n , b n := u n - iw n , 

In ■ — "n ~r ^"nj 9n ■ — "n ^^m 



nez. 



The inhomogeneous system (|13p transforms to the equivalent form 



m£U+UU- u n,m(\*- ~^~ P) a m + P°m) 1'fln — Jni 



-e(a n +i — 2a n + a n _i) + a n — J2meu + uU- ^ 

-e(b n+ i - 2b n + b n -x) + b n - ^Z meU uu _ 5 n ^ m (pa m + (1 + p)b m ) + Q.b n = g n , 



(14) 



which can be rewritten in the operator form 



L 



a 




a 




f 




b 


-0 


b 


— 


. -S _ 


, L = 



-eA + I-(l + p)V -pV 

P V eA-I+(l + p)V, 



(15) 



where A : l 2 ( 



Z 2 (Z) is the discrete Laplacian operator 



(A«) n := u n+1 - 2u n + n n _i, n E 



and V : Z 2 (Z) — > l 2 (jL) is the associated compact potential 



(Vu) n = E 5 n . 

m&U+UU- 



n E 



Let Rq(\) : Z 2 (Z) — > Z 2 (Z) be a free resolvent of the discrete Schrodinger operator —A for 
A ^ o~(— A) = [0,4]. The free resolvent was studied recently by Komech, Kopylova, & Kunze |14j . 
The free resolvent operator can be expressed in the Green function form 



Vf G / 2 (Z) : (i2o(A)f) T 



1 



£ 



-iz(\)\n—m\ 



2isinz(A) 

where z{\) is a unique solution of the transcendental equation for A ^ [0,4] 

2-2cosz(A) = A, Rez(A) G [-7r,vr), Imz(A) < 0. 



J mi 



(16) 



(17) 



The limiting absorption principle (see, e.g., Pelinovsky & Stefanov |19| ) states that a bounded 
operator R (X) : l 2 (Z) ->■ / 2 (Z) for A £ [0,4] admits the limits 



Rq(oj) = lim i? (w ± ifi) : g (Z) -4- lL ff (Z) , cr > 



/40 



for any fixed w G (0, 4). 

The limiting free resolvent operators R Q (w) can also be expressed in the Green function form 



Vf G l\Z) : (<(u;)f) n 



V^ e ^ i0 ±( uj )\ n ~ m \ f 



(18) 



2i sin 0-t (a;) 

31 v ; mez 

where 0±(u;) = ±6(uj) and 0(w) is a unique solution of the transcendental equation for uj G [0,4] 

2-2cos6(uj)=uj, Re6(uj) G [— 7r,0], Im0(w) = 0. (19) 

The limiting operators i? (uj) : / X (Z) — )■ Z°°(Z) are bounded for any fixed uj G (0,4) but diverge as 
U) 4 and w|4 These divergences follow from the Puiseux expansion 



vf g i\{i) : (n£(u>)f) n = ^TT E /- - ^ E l n - m l/m + (^M f )n. 



(20) 



where 



3C > : [|A±(w)f||,oo < C|0 ± (w)|||f|| / i. 



Divergences of R (uj) at the end points uj = and w = 4 indicate resonances, which may result 
in the bifurcation of new eigenvalues from the continuous spectrum on [0, 4] either for A < or 
A > 4, when —A is perturbed by a small potential in / 2 (Z). 

Let us denote the solution of the inhomogeneous system (|15p by 



a 
b 



Rl{V) 



g 



R L (ty 



Rn(n) i2i 2 (fi) 

«2i(n) i? 22 (^) 



(21) 



The following theorem represents the main result of this section. This theorem is valid for the 
simply connected sets U + U U- , which are defined by the following definition. 



Definition 2 We say that the set U+ U U- is simply connected if no elements in Z\{£/+ U £/_} are 
located between elements in U+ U U- . 

Theorem 1 Fix U+, U- C Z suc/i i/tai U+C\U- = 0, N := \U+\ + \U-\ < oo, and U+ U L7_ is 
simply connected. There exist small eo > and 5 > suc/t t/iai /or any fixed e £ (0, eo) the resolvent 
operator 

R L (Q) : l 2 (Z) x l 2 (Z) -> Z 2 (Z) x Z 2 (Z) 

is bounded for any Q ^ -B<s(0) U [1, 1 + 4e] U [— 1 — 4e, —1]. Moreover, Rl(£1) has exactly 2N poles 
(counting multiplicities) inside Bg(0) and admits the limits 

RJ;(P) :=HmR L (n±in) 



such that for any £7 G [1,1+ 4e] U [—1 — 4e, —1] and any e G (0, eo), there is C > such that 

>L '-'-yiUjxiJ-^xZ 00 



|i? L (ri)||^l v ;l^;oo v /oo <Ce 



Remark 3 The other way to formulate the main theorem is to say that the end points of the 
continuous spectrum a c {L) = [1,1 + 4e] U [— 1 — 4e, —1] are not resonances and no eigenvalues of the 
linear operator L may exist outside a small disk Bs(0) C C. The 2N eigenvalues inside the small 
disk Bg(0) are characterized in Proposition^ 

Solving the linear system (I14p with the Green function (I16p . we obtain the exact solution for 
any n £ Z 



2iesinz(A^ 
1 



E 






-iz(Xj^)\n—m\ 



fm + E 



radU+UU- 



-iz{\ + )\n- m \^ l+p)am+phm) " 



bn = ^hnr^ ( Emez e'^-^-^Qm + £ meC/+uC7 _ e -<*(A-)l"-"» \(pa m + (1 +p)b m ) ) , 



(22) 



2ies'mz(\-) 

where the map CsA^zfCis defined by the transcendental equation (|17p and 

±0-1 



e 

The solution is closed if the set {(a n , o n )}ne(7+u!7- is found from the linear system of finitely many 
equations for any n G U+ U U- 



2iesinz(A + )a„ — £ 7 
2ie sin z(\^)b n -J2 m eu, 



mdU+VJU. 



UU- 



e -iz{\ + )\n-m\((l +p)am+phm) = ^ mgze - 4 ,(A + )|n- m | /m) 



-iz(A_)|n— m\ 



(pa m + (1 + p)b m ) = Em&z e 



-iz(\—)\n— m\ 



fjn 



(23) 



Let us order lattice sites n G U+ U U- such that the first site is placed at n = 0, the second 
site is placed at mi, the third site is placed at mi + m2, and so on, the last site is placed at 
mi + m i + • • • + mjv-l, where iV = |C/+| + |C/_| and all ?rij > 0. If [/+ U f7_ is a simply-connected 
set, then all rrij = 1. 

Let Q(qi, g2, • • • , ojv-i) be the matrix in C defined by 



Q(Qi,Q2,--- ,Qn-i) 



1 


91 


9192 


• • • 9192 • 


• • 9JV-1 


91 


1 


(12 


■ ■ ■ 9293 • 


• • ?iV-l 


9l92 


92 


1 


... gg . . 


•gjv-i 



9192---9AT-1 92---9Af-i 93---9JV-1 



(24) 
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Figure 1: Schematic display of various domains in the 0-plane. 



Let qf = e m 3 z ( x ±) an d Q ± (Q,e) := Q(qf,qf,--- ,Q^r i)- The coefficient matrix of the linear 
system (|23p is given by 



>Qn-1' 



A(Q,e) r- 



2iesinz(\ + )I - (1 + p)Q+(0, e) ~pQ + (^, e) 

-pQ-(n.e) 2»esin*(A_)I-(l+p)Q~(n,e) 



(25) 



where I is an identity matrix in C . 

We split the proof of Theorem Q] into three subsections, where solutions of system (|22p and ([23 
are studied for different values of f2. 



3.1 Resolvent outside the continuous spectrum 

We consider the resolvent operator Rl(Q) for a fixed small e G (0, eo). The following lemma shows 
that Rl(£1) is a bounded operator from l 2 (7L) x 1 2 {7L) to l 2 (7L) x l 2 (Z) for all SI G C except three 
disks of small radii centered at {0, 1, —1}. 

Lemma 1 There are eo > and S,6± > such that for any e G (0,eo), the resolvent operator 
R L (Q) : l 2 {%) x l 2 (1) -> Z 2 (Z) x Z 2 (Z) is founded /or all tt G C\{B,$(0) U Sj + (1) U 5 5 _(-l)}. 
Moreover, Rl(Q) has exactly 2N poles (counting multiplicities) inside B$(0). 

Proof. From the property of the free resolvent operator Rq(X), we know that the Green function 
in the representation (|22p is bounded and exponentially decaying as |n| — > oo for any SI such that 
A± i [0,4]. This gives Q, £ a c (L) = [1,1 + 4e] U [-1 - 4e, -1]. Therefore, R L (n) is bounded map 
from l 2 (Z) x l 2 (Z) to l 2 (Z) x / 2 (Z) for any ^ a c (L) if and only if the system of linear equations ([23]) 
is uniquely solvable. We shall now consider the invertibility of the coefficient matrix A(£l, e) of the 
linear system ([23|) in various domains in the $7-plane for small e > 0. Figure Q] shows schematically 
the location of these domains on the O-plane. 

Fix e)o G (0, 1). Let O belong to the vertical strip 

S 5o :={QeC: Re(n)e[-S ,S }}. 

Then z(X±) = —in± are uniquely determined from the equation 

l + O 



K± 



+ e 



-K± 



-, Re(/s±) > 0, Im(K±) G [— tt, it), 



,4(0,0) 



(26) 



which admits the asymptotic expansion 

e 1=f^ 

and 

e~ K± = — e — + (D{e 2 ) as e ->• 0. 
1 =F^ 

Therefore, both esinh(K±) and Q^(Q,e) are analytic in e near e = and 

2zesinz(A ± ) = 2esinh(K ± ) = l=Fil + 2e + C(e 2 ) as e ->• 

and 

Q ± (n,e) = I + 0(e) as e -»• 0. 

It becomes now clear that A(Q, e) is analytic in O G S$ and e G (— eo, eo) with the limit 

-(p + n)i -p/ 
_ p i -(p-o)i 

Matrix ,4(0, 0) G C 2Nx2N is singular only for O = 0. Thanks to analyticity of ,4(0, e), the 
determinant D(0,e) = det,4(0,e) is also analytic in these variables and 

D(n,e) = (-n 2 ) N + 0(e) as e ->■ 0. 

Therefore, there exist 2N zeros of .D(0,e) for small e G (0, eo) in a small disk -B«5(0) with 5 = 
0(e 1 ' 2N ). By Cramer's rule, these zeros of D(0,e) give poles of Rl(Q). 
Fix 5+ G (0, 1) and + G (f ,7r). We now consider O in the domain 

S 6+ :={n = l + re ie , r > 6+, fleH + ,0 + )}. 

In this domain, we have the same presentation for #(A_) = — in- but a different presentation for 
z(A + ) = —ik + — ir. Now k_|_ is uniquely determined from the equation 

e K + + e~ K + + 2 = ^—^ = -e ie , Re(n + ) > 0, Im(/s + ) G [0, 2tt), 
e e 

which admits the asymptotic expansions 

e K + = ?—± - 2 - -^-j + 0(e 2 ) as e^0 

and 

2ie sin z(X+) = 1 - + 2e + C(e 2 ) as e ->• 0. 

Since Re(K+) — > oo as e — >■ 0, .4(0, 0) is the same as matrix (|26p and it is invertible for O G S < 5 + . 
Similar arguments can be developed for 

S s „ := in = -1 + re ie , r > 5-, 9 G (0_, 2ir - 9-)\ , 

where 5- G (0, 1) and 9- G (0, ^). Because there are choices of 5o,5± > such that 

S,5 U S s+ U S*_ = C\{B S+ (1) U S 5 _ (-1)}, 
we obtain the assertion of the lemma. ■ 

Remark 4 The proof of LemmaUl implies that poles of Rl(Q) may have size |0| = 0(e 1 ' ). The 
results of the perturbation expansions (see )11^ for details) imply that the eigenvalues bifurcating 
from in the full spectral problem (£5|) have size 0(e 1 ' 2 ). Moreover, the same perturbation expansion 
technique can be applied to show that eigenvalues of the truncated spectral problem M3\) have the 
same size 0(e 1 ' 2 ). 

9 



3.2 Resolvent inside the continuous spectrum 

We shall now consider the resolvent operator Rl(Q) inside the continuous spectrum 

<7 c (L):=[l,l + 4e]U[-l-4e,-l]. 

Thanks to the symmetry of system (|22j) - (123p in fi, we can consider only one branch of the continuous 
spectrum [1, 1 + 4e]. Therefore, we set = 1 + eui with co £ [0,4] and define 

z(A+) = z(uj) = 9 and z(A_) = z{— 2e~ — u>) = —in. 

It follows from (|17p and (jT9|) that 9 £ [— it, 0] and k > are uniquely defined from equations 

2-2cos(0) =u, 2e(cosh(K) - 1) = 2 + eu, we [0,4]. (27) 

The choice of 9 £ [— n, 0] corresponds to the limiting operator R^(ui) of the free resolvent. Since 
Rq(u) : IffCZ) —7- /^(Z) is well defined for u £ (0,4) and a > tj, -^(H~ ew ) is a bounded map from 
/ 2 (Z) x /^(Z) to /^ CT (Z) x /^(Z) for any uj £ (0,4) and ct > \ if and only if there exists a unique 
solution of the linear system (f23|) . On the other hand, the free resolvent is singular in the limits 
U) \. and w|4 and, therefore, we need to be careful in solving system ([22l) - (123p in this limit. 

The main result of this section is given by the following theorem. 

Theorem 2 Let m\ = rri2 = ■ ■ ■ m^^i = 1. There exists eo > such that for any to £ [0,4] and 
any e £ (0, eo), there exist C > such that 

\\R+(1 + eu)\\ llxll _ laoxlao < Ce-\ (28) 

where the upper sign indicates that ui is parameterized by uj = 2 — 2cos(#) for 9 £ [— it, 0]. 

To prove Theorem we analyze solutions of system (123|) for uj £ [0,4]. Let us rewrite explicitly 

q+ = e- im i e and qj = e"" 1 ^, j £ {1, 2, ..., N - 1}. 

The coefficient matrix (|25p for (7 = 1 + eoj with uj £ [0,4] is rewritten in the form 



A(9,e) 



2ie sin(9)I - (1 + p) Af (0) - P M(9) 

-pN(n) 2e sinh(K)/ - (1 + p)N(k) 



(29) 



where M{9) = Q(qf ,q£ , • • • ><zjv-i) anc ^ N(k) = Q(q^ , q^ , • • • i?^]^). Note that and M(#) are 
e-independent, whereas N(k) depends on e via k. The linear system (123|) is now expressed in the 
matrix form 

A(0,e)c = fc(0,e), (30) 

where components of c £ C and /i £ C are given by 

1 f V p -i0|ri-m| f 



&n J neU+UU- I S me z e "" m fe ) n&U+YJU- 

Thanks to the asymptotic expansion 

2 e n 

e K = - + 2 + w-- + C(e 2 ) as e -> 0, 

we have 

2esinh(K) = 2 + (2 + w)e + 0(e 2 ) as e -»• 0. 

Both -4(0, e) and h(9, e) are analytic in 9 £ [—it, 0] and e G (— eo, £o)- The following lemma establishes 
the invertibility condition for matrix A{9, e). 
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Lemma 2 For any e G (0, eo), matrix A(9,e) has a zero eigenvalue of geometric and algebraic 
multiplicities N — 1 for 9 = —it and 9 = 0. If mi = 777,2 = ••• = m^-i = 1, matrix A(6,e) is 
invertible for any 9 G (— it, 0). 

Proof. We use the fact that matrix A(9,e) is analytic in e for small e G (— eo,£o)- Therefore, it 
remains invertible if A(9, 0) is invertible. To consider the limit e — > 0, we note that k — > 00 and 
N(k) — > I as e — > 0, so we have 



A(0,O) 



{l+ P )M{6) -pM(9) 
—pi (1 — p)I 



For any p € N, matrix A(9, 0) is invertible if and only if matrix M(9) is invertible. Let us then 
compute 

DN(qi,q2,--- ,QN-i) ■= detQ(q 1 ,q 2 , ■ ■ ■ ,<?jv-i). 

and Dj\i(qi,q2, ■ ■ ■ , qN-i) is a quadratic polynomial of q\. 



We note that Djy(±l, q 2 , ■ ■ ■ , qN-i 
Therefore, 

DN{qi,Q2, • • • , ^at-i) = (1 - qf)D N (0, q 2 ,--- , qN-i) = (1 - ql)D N _i(q 2 , ■ 

Continuing the expansion recursively, we obtain the exact formula 

D N (q 1 ,q 2 ,--- ,q N -i) = {1 - qf){l - <?f) • • • (1 - g^_i), 



.qn-i, 



(31) 



from which we conclude that Q(qi,q 2 ,- ■ ■ ,qN-i) is invertible if and only if all g^ 7^ ±1. This 



implies that M(9) is invertible if and only if all e 



7^ ±1, which is satisfied if all nij = 1 and 



9 G (— vr,0). The second assertion of the lemma is proved: for any e G [0,eo), matrix A(9,e) is 
invertible for 9 G (— 7r, 0) if all 771^ = 1. 

The first assertion of the lemma tells us that for any e G (0,eo), matrices A + (e) := ^4(0, e) and 
A_{e) := A(—tt, e) have a zero eigenvalue of geometric and algebraic multiplicities N — 1. We write 
.A-i- (e) explicitly in the form 



A±(e) 



[l+p)M± -pM± 

-pN(K±) 2e sinh(K±)I - (1 + p)N(k±) 



where k± > are uniquely defined by 

2e(cosh(K + ) — 1) 
whereas matrices M± are given by 



M 4 



2, 2e(cosh(K_)- 1) = 2 + 4e, 



111- 


• 1 


111- 


• 1 


111- 


• 1 



1 1 1 



and 



M_ 



1 


(-1)" 11 


C_l~) m l+ m 2 


... (-y\m 1 +rn2 + ...+m N 


(-l) ml 


1 


(_l)m 2 


t_-\\m2 + ...+m N --L 


_2) m l+ m 2 


(-1)™ 2 


1 


C_l) m 3 + ---+"iiV-i 



( , _]_ s ) m l+ m 2+---+"'-Ar-i (_-|\m2+...+m i v_i /'_^\m3+...+m iv _ 
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It is clear that Null(M+) and Null(M_) are (N - l)-dimensional. 

The first N rows of A^(e) are identical to the first row, whereas the last N rows of A + {e) are 
linearly independent at e = and, by continuity, for small e G [0, eo). Therefore, Null(A + (e)) is 
[N — l)-dimensional for any e G [0, eo)- Similarly, the second, third, and JV-th rows of j4_(e) are 
identical to the first row multiplied by (-l) m \ (-l)™i+™2, and (_1)™i+™2+-+™jv-i respectively. 
The last N rows of ^4_(e) are linearly independent for small e > 0. Therefore, Null(A_(e)) is 
(N — l)-dimensional for any e G [0, eo). 

It remains to prove that the zero eigenvalue of A± (e) is not degenerate (has equal geometric and 
algebraic multiplicity) for e G (0, eo). It is clear from the explicit form of A±(0) and M± that 



t/GNull(4(0)) ^ 



(1 — p)w 
pw 



w G Null(M±). (32) 



To construct a generalized kernel, we consider the inhomogeneous equation 

A±(0)u = u, ueNull(At(0)). 
Then, we obtain for w G Null(iW±), 



u 



[1 — p)w — w 
pw 



M±w = (jp— l)w. 



lip 7^ 1, then no w G C^ exists because M± is symmetric. Therefore, for p^ 1, the zero eigenvalue 
has equal geometric and algebraic multiplicity for the matrix A±(0) and, by continuity, for the 
matrix A±(e) for e G [0, eo). 

The case p = 1 needs a separate consideration since w = and the zero eigenvalue of A± (0) has 
geometric multiplicity N — 1 and algebraic multiplicity 2N — 2. This case is considered in Appendix 
A, where we show that the degeneracy is broken for any e 7^ 0, so that A±(e) in the case p = 1 still 
has a zero eigenvalue of equal geometric and algebraic multiplicity A^ — 1 for any e G (0, eo). ■ 

Because the coefficient matrix A(9, e) is singular at 9 = and 9 = —tv, we shall consider the 
limiting behavior of solutions of the linear system (130j) near these points. The following abstract 
lemma gives the sufficient condition that the unique solution c of the linear system (|30p for small 
9 7^ and fixed e G (0, eo) remains bounded in the limit 9 — > 0. Because e is fixed, we can drop this 
parameter from the notations of the lemma. 

Lemma 3 Assume that A(9) G C 2Nx2N and h{9) G C 2N are analytic in 9 G (-9 ,9 ) for 9 > 
and consider solutions of 

A(9)c = h(9), ceC 2N . 

Assume that A{9) is invertible for 9^0 and singular for 9 = and that the zero eigenvalue of A(0) 
has equal geometric and algebraic multiplicity n < 2N . A unique solution c for 9 7^ is bounded as 
9^0 if 

h(0) ±Null(A*(0)) and Null^O)^^))) = {0}- (33) 

Remark 5 We denote the Hermite conjugate of a matrix Aq G C x by Aq = Aq . Let 

Null(.Ao) = span{«i, ...,u n } and Null(-Ao) = spanjui, ...,v n }, (34) 

where {u\, ...,u n } and {v±, ...,v n } are mutually orthogonal bases, so that 

(ui,Vj) C 2N = 5i t j for all 1 < i,j < n. (35) 

12 



The restriction of matrix A\ £ (£ 2Nx2N on Null(Ao) denoted by -AiInuII(Ao) can ^ e expressed by the 
matrix P G <£ nxn w ah elements 

Pij = (vj, Aim) C 2N for all 1 < i,j < n. (36) 

Proof. The proof of the lemma is achieved with the method of Lyapunov-Schmidt reductions. 
Using analyticity of A(9) and h(9), let us expand 

A{9) = A + OAi + 9 2 A{9), h{9) = h + 0/ti + 9 2 h(9), 

where A = ,4(0) , A\ = A'(0), h = h(0), hi = /i'(0), and A(9) and h(9) are bounded as 9 ->■ 0. 
Given the basis for Nu11(j4o) in ([34"[) . we consider the orthogonal decomposition of the solution 

n 

c = Y^ ajUj + b, b± Null(Ao). (37) 

i=i 

The linear system becomes 

n 

(A) + 0A 1 + 9 2 A{9))b + 9j2 a j ( A i + 0A{9)) Uj = h + 9h x + 9 2 h{9). (38) 

i=i 

Projections of system (I38p to the basis for Nu11(^4q) in (I34p give n equations 

n 

Y, (Pij + 0Pij{0)) dj + (vi, (A ± + 9A(9))b) C 2N = (v i} hi + 9h{9)) C 2 N , 1 < i < n, (39) 

3=1 

where Pij is given in (|36p . Pij{9) = (vi,A(9)uj) C 2N is bounded as 9 — > 0, and we have used the 
condition h _L NuII(Aq). 

Let Q : C 2N ->■ Ran(A ) C C 2N and Q* : C 2N ->■ Ran(A^) C C 2N be the projection operators. 
Recall that Ran(A ) 1 Null(^) and Ran(^) _L Null (4)). Projection of system ([38]) to Ran(A ) 
gives an equation for b 

n 

Q(A + 9Ai + 9 2 A{9))Q*b = Q{h + 9h x + 9 2 h{9)) - ^ a J Q(A 1 + 9A{9))u r (40) 

i=i 

Because QAqQ* is invertible, there is a unique map C™ 9 (ai,...,o n ) i— >■ 6 € Ran(Ag) for any 
€ (—#0) $o) such that b is a solution of system ([40 p and for any 9 G (— 9q, 9q), there is C > such 
that 

||b - QM^Q/iolbiv < C9. (41) 



Since Nu11(^4i|nu11(A )) = {0}, matrix P is invertible. For any b from solution of system (|40p 
satisfying bound (|4"T|) . there exists a unique solution of system ([39]) for (oi,...,a n ) for any € 
{—9q,9q) such that 

3C>0: ||a-P- 1 (/-Q)(/ i i-AiQM 1 Q/i )||c™ < C0. (42) 

For any 0^0, the solution of system A(9)c = h{9) is unique. Therefore, the unique solution 
obtained from the decomposition ([37]) for any 9 £ (—9q,9q) is equivalent to the unique solution of 
system ,4(0)c = /t(0) for 0/0. ■ 
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We shall check that the conditions (|33p of Lemma [3] are satisfied for our particular matrix A(9, e) 
and the right-hand-side vector h(9, e) for both end points 9 = and 9 = —n. 

Lemma 4 Let h + (e) := /i(0,e) and /i_(e) := h(—7r,e). For any e G (0, eo), it is true that 

h±(e)±mil(A* ± (e)) and Null(cM±(e)| Null(y4±(e)) ) = {0}. (43) 



Proof. It is sufficient to develop the proof for = 0. The proof for 9 = —n is similar. 

Recall that the first N rows of A(0, e) are identical to the first row. Since components of h(0, e) 
are given by 



E 



n\n-m\ ( 

Urn J neU+UU- 



the first N entries of h(0, e) are also identical so that h(0, e) G Ran(^4(0, e)) _L Nu11(j4*(0, e)) for any 
e G (0, eo). Therefore, the first condition ([43 1) is satisfied. 
Next, we compute -Ai(e) = d$A(9, e)|e = o- We know that 



2e(cosh(«) - 1) = 2 + e(2 - 2 cos(0)) 



c?k sin(0) 



ci# sinh(K) 



therefore, 



where 



M*) 



2eI+{l+p)R pR 




(44) 



R 



77T-1 


mi + m 2 


• • • mi + m 2 H h mjv-i 





m 2 


• • • m 2 + m 3 -| h mjv_i 


m 2 





m 3 -\ h m N ^i 




mi 

m-i + m 2 



mi + m 2 H h mAr_i m 2 -I h mAr_i m 3 -\ + m N ^i 



Let -P(e) be the matrix in C' l)x(« i) wrr j cri represents the restriction Ai(e) 



Null(A (e))- 



Ex- 



istence of a G Null(P(e)) C C is equivalent to existence of u G Null(vlo(e)) C C such that 

Ai(e)u G Ran(A (e)) _L Null(>^(e)). In other words, we need to find u G Null(A)(e)) such that the 
first N entries of A\(e)u are identical (the other N entries of A\(e)u are zeros). 

By continuity in e G [0, eo), the second condition (143 P is satisfied if it is satisfied for e = 0. 
Therefore, it is sufficient to check the existence of u G Nu11(j4o(0)) such that the first N entries of 
Ai(0)u are identical. 

It follows from relations (j32p and (|44p that existence of u G Null(.Ao(0)) such that the first N 
entries of Ai(0)u are identical is equivalent to the existence of w G Null(M+) C C such that all 
entries of Rw are identical. 

If w= [w 1 ,w 2 ,-..,w N ] T G Null(M + ), then 



'll'i 



Condition (Rw)i = (Rw)?, gives 



+ W 2 + ... +WN = 0. 



(45) 



m\(w2 + ••• + wn) = m\w\. 
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Constraint (|45p implies that if mi 7^ 0. then w\ = and W2 + • • • + Wjv = 0. Continuing by 
induction for condition (Rw)j = (Rw)j+i, where j £ {1, 2, ...,N — 1}, we obtain that if rrij ^ 0, 
then Wj = for all j £ {1, 2, ..., N — 1}. In view of constraint (|45|) . we have wn = that is 
w = £ C . As a result, we have proved that Nu11(-Ai(0)|nu11(A (o))) = {0}- By continuity in 
e £ [0, eo), Null(^4i(e)|Null( J 4 (e))) = {0} f° r sman e 7^ 0, which gives the second condition f|43[) for 
= 0. ■ 

Remark 6 Lemma^is proved without assuming that all rrij = 1. 

Proof of Theorem [2j By Lemma [H assumptions of Lemma [3] are satisfied and the unique solution 
of system (I30p for £ (— 71", 0) is continued to the unique bounded limit cq = lim^o c - From the 
first iV equations of system (I23p , we infer that 



6 = 0: Yl ((I + P)a m + pb m ) = ~ ^ fm- 

m&U+VJU- meZ 

As a result, the simple pole singularity at 6 = (z(A+) = 0) in the Green's function representation 
(|22p with the Puiseux expansion (|20p is canceled. Similarly, the simple pole singularity at 9 = —ir 
is cancelled. On the other hand, the representation ([22j) contains e in the denominator, which does 
not cancel out generally. As a result, Lemma [2] for all rrij = 1 and Lemma U] give that for any 
to £ [0,4] and any e £ (0, eo), there exists C > such that 

||a||/oo < Ce _1 . 

This gives bound (J28l) and hence Theorem [2j □ 

3.3 Matching conditions for the resolvent operator 

To complete the proof of Theorem [TJ we need to prove that no singularities of linear system (J23H 
are located inside the disks Bs + (1) and Bg_(—1) for e-independent 6± > 0. It is again sufficient to 
consider the disk Bg + (1) because of the symmetry in the fi-plane. 

The free resolvent operator Rq(X) : ^(Z) ~~ ^ ^— <r(^) with a > \ is extended meromorphically 
in variable (9(A) for A £ C+\[0,4] with simple poles at 6 = (A = 0) and 9 = —k (A = 4). By 
Theorem El the resolvent operator R+(l + eu) : *}(Z) x Z*(Z) -)• Z°°(Z) x /°°(Z) is bounded for 
a; £ [0,4] and the pole singularities are canceled. As a result, the resolvent operator -R^(l + eA) 
can be extended as a bounded operator from Z 2 (Z) x Z 2 (Z) to /^(Z) x /^(Z) with c > 5 for any 
A £ C + \[0, 4]. We need to show that no singularities of the resolvent operator i?i(l + eA) exist in 
the upper semi-annulus 

D s+ = {\eC+: 7+ <\X\<5 + er 1 }cB 5+ (l), 

where 7+ > 4 and 5+ £ (0,1). A similar analysis can also be used to show that the resolvent 
operator i?^(l + eA) can be extended as a bounded operator in the lower semi-disk in Bg (1). 

Lemma 5 For any e £ (0, eo) and all A £ Dg + , the resolvent operator i?i(l + eA) is a bounded 
operator from Z 2 (Z) x Z 2 (Z) to Z 2 (Z) x Z 2 (Z). 

Proof. Since the continuous spectrum does not touch boundaries of D§ , the statement is true if 
and only if there exists a unique solution of linear system (|23p . 
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Let us denote z(A-f) = z(X) and z(X—) = —in(X), where z(X) is found from the transcendental 
equation (JTTj) and k(A) with Re(ft(A)) > admits the asymptotic expansion for A G D$, 



,/c(A) 



2 + eA 



+ 2 



2 + eA 



+ C(e^) as e^O. 



As earlier, we denote qt = e im i z ( x ) and q- = e m o K ( x ) for j G {1, 2, ..., iV - 1}. 
We write the coefficient matrix (1251) for il = 1 + eA in the form 



A(X,e) 



-pM(\) 



-e v / A(A-4)I-(l+p)M(A) 

-piV(K) ^(2 + eA) 2 + 4e(2 + e\)I - (1 + p)JV(«) 



(46) 



where M(A) = Q(qf , q^ , • • • , <Zjv_i)> N{n(X)) = Q{q\ , <? 2 > ' ' ' > 9jv-i)> an( ^ * ne appropriate branches 
of sinz(A) and sinh(«(A)) are chosen in the domain -D,5 + . 
Let |A| = e>(e~ r ) as e -)• for r G [0, 1). Then, we have 



A(X,e) 



-(l+p)M(A) -pM(A) 
-pi (l-p)I 



as e — > 0, 



(47) 



where M(A) — > I as e — > if r G (0,1) and M (A) ^> I as e — >■ if r = 0. The limiting matrix (I47p is 
not singular if 7+ > 4. Hence A(A, e) is not singular for small e > if |A| = 0(e _r ) with r G [0, 1). 
Let |A| = C(e~ r ) as e -)• for r G (0, 1]. Then, we have 



A(X,e) 



[l + eX + p)I -pi 

-pi (l + eX-p)I 



as e 



0. 



Again, the limiting matrix is not singular if eA 7^ —1 (that is <5+ < 1) and hence A(X, e) is not 
singular for small e > if |A| = 0(e~ r ) with r £ (0, 1]. 

Since the above asymptotic scaling overlap at r G (0,1), the matrix A(X,e) is not singular in 
the domain D& + for small e > 0. ■ 

Theorem Q] is proved with Lemma Q] Theorem [2] and Lemma [5] 



4 Perturbation arguments for the full resolvent 

Let us now consider the full spectral problem ([5]). Thanks to Proposition Q] and expansion (J7J), we 
can represent (ftn by 

meU + l)U- 

where {Xm}m£U+uU- 1S a se t of numerical coefficients and {W n } n ^z £ l 2 (Z) is a new potential such 
that ||W||j2 = O(l) as e -)■ 0. 

In variables {(a n ,fr n )} ne z> the resolvent problem can be rewritten in the operator form 



where 



L 



(L + e 2 H>) 



-eA + I -(l+p)V -pV 

pV eA-I + (l+p)V, 



a 
b 


-n 


a 
b 


= 


f 



(49) 



W 



'l+p)W -pW 
pW {l+p)W, 
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and V is the associated compact potential such that 

(Vu) n = 



meU+ut/_ 



Let us denote the solution of the inhomogeneous system (|49p by 



a 

b 



R(n) 



f 
-g 



(50) 



where i?(f2) is the resolvent operator of the full spectral problem ([8]). The following theorem 
represents the main result of our paper. 

Theorem 3 Fix U + , C/_CZ such that U + C\U- = 0, N := \U+\ + \U-\ < oo 7 and U+UU- is simply 
connected. For any integer p > 2, there are eo > and 5 > such that for any fixed e £ (0, eo) the 
resolvent operator 

R(n) : / 2 (Z) x Z 2 (Z) -> / 2 (Z) x l 2 (Z) 

is bounded for any Q ^ B$(0) U [1, 1 + 4e] U [—1 — 4e, —1]. Moreover, R(il) has exactly 2N poles 
(counting multiplicities) inside B$(0) and admits the limits 

R ± (n) :=hmR(n±in) 

/aJ.0 

such that for any f2 £ [1,1 + 4e] U [—1 — 4e, —1] and any e £ (0, eo), there is C > suc/i £/iat 



ff*(n)[| 



Zjx2j-^°°xi c 



<Ce 



-i 



Proof. Let R^(Q) be the resolvent operator for the inverse operator (L — fi/) -1 associated 
with the compactly supported potential V. We shall prove that Theorem [1] remains valid for the 
resolvent operator R^(JT). Assuming it, the rest of the proof of Theorem [3] relies on the perturbation 
arguments and the resolvent identities 

R(n) = R L (n)(i + e 2 wR~ L {£i))- 1 = (/+ e 2 R L (n)wy 1 R L (n). 

Indeed, outside the continuous spectrum located at 

a c (L + e 2 W) = a c {L) = a c {L) = [-1 - 4e, -1] U [1, 1 + 4e], 

the resolvent operator R^(Q) is only singular inside the disk Bg (0), where perturbation theory 
of isolated eigenvalues apply. Inside the continuous spectrum, R^(Q) is extended as a bounded 
operator from |J(Z) x Z X (Z) to Z°°(Z) x /°°(Z) such that for any O £ [1, 1 + 4e] and any e £ (0, e ), 
there is C > such that 



3c>0: \\Rf(n)\\ 



llxll^l^xl 00 



<Ce 



(51) 



Since W is a bounded (f2,e)-independent operator from Z°°(Z) x Z°°(Z) to Zj(Z) x /J(Z) (note here 
that £ 1 2 , 2 (7j), see Remark [T]), bound (J51D implies that 

3C>0: \\e 2 WR L (n)\\ llxll ^i xll <Ce, 
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so that (I + e 2 WR~ L (Q)) is an invertible bounded operator from l\ (Z) x /{(Z) to l\ (Z) x Z{(Z) for 
small e > 0. 

We only need to extend Theorem Q] to the resolvent operator R^(Q). The Green's function 
representation (122|) and the linear system (J23l) are now written with the factor (1 + eXm) in the sum 
over m £ U+ U U-. This implies that the coefficient matrix A(Q, e) is now written as 



A(Q,e) :-- 



2iesin z(X+)I - (1 + p)Q + (tt, e)(J + eD) -p<2 + (^, e)(J + eD) 

-pQ~(n,e)(I + eD) 2iesinz(A_)I - (1 + p)Q~(tt,e)(I + eD) 



where D is a diagonal matrix of elements {Xm}meU+uU-- If P > 2, Lemmas [I] [2] [H and [5] remain 
valid as these lemmas were proved from the limit e = (perturbation theory of Appendix A is only 
required for p = 1), where A(£l, 0) = A(£l, 0). Therefore, Theorem [1] holds for the resolvent operator 
R L (n)iip>2. ■ 

Corollary 1 The result of Theorem holds for p = 1 if N = 1 . 

Proof. If A^ = 1 (which is the case of fundamental discrete soliton). the 2x2 coefficient matrix 

A(Q,e) ■ 



2iesinz(A+) - (l+p)(l + exo) -p(l + e%o) 

-p(l + exo) 2iesinz(A_) - (1 +p)(l + exo) 



is only singular in 5,5(0) for small e > 0, where a double pole of R^(£l) and R(£l) resides. ■ 

Unfortunately, in the cubic case p = 1, we can not generally extend the result of Theorem [3] to 
multi-site discrete solitons with N > 2 because the perturbation theory for A(£l, e) near the end 
points of the continuous spectrum fi = ±1 and fi = ±(1 + 4e) draws no conclusion in a general 
case. For instance, reworking the perturbative arguments of Appendix A, we obtain the necessary 
condition for Null(A-|-(e)) 2 > Null( J 4 ± (e)) in the form 

e(2/ - J - 2D)w + 0(e 2 ) ±w£ Null(M + ), 

where / is the identity matrix in K . , J is the two-diagonal matrix (f55|) from Appendix A, and D 
is a diagonal matrix of {Xm}meU+uU- ■ Because (21 — J — 2D) is no longer positive definite, the 
degenerate cases with Null(^4-|-(e)) 2 > Null(A-|-(e)) are possible. 

To illustrate this possibility, we set N = 3 and consider three distinct simply-connected discrete 
solitons associated with the sets 

(a) 17+ = {0,1, 2}; (b) U+ = {0, 1}, U-{2}; (c) U+ = {0,2}, 17_{1}. 

Computations of the power expansions ([7]) give 



(a) Xr, 




(b) X, 




( c ) Xr. 




As a result, matrix C — 21 — J — 2D is obtained in the form 



(a)C 



-1 
-1 2 
-1 



(b)C 



-1 
-2 -1 
-1 -4 



(c)C 



4 


-1 





1 


-6 


-1 


) 


-1 


-4 



We have 



Null(M_|_) = span{wi,W2}, w\ 



V2 



W-2 



1 



from 



which we compute the matrix of projections P^ = (Cwi,Wj)t£3 in the form 



(a)P 








3 J 



, (b)P 



2 



2 

5 

3 



(C)P 







The projection matrices in cases (a) and (b) are singular. In order to show that Null(74-|-(e)) 2 = 
Null(j4±(e)) for e G (0, eo), we need to extend perturbation arguments of Appendix A to the order 
0(e 2 ). Although it is quite possible that the non-degeneracy condition Null(j4-|-(e)) 2 = Null(.A±(e)) 
is still satisfied for simply-connected multi-site discrete solitons for p = 1, we do not include com- 
putations of the higher-order perturbation theory in this paper. 

5 Case study for a non-simply-connected two-site soliton 

We explain now why the resolvent operator associated with non-simply-connected multi-site discrete 
solitons have singularities near the anti-continuum limit. These singularities appear in Lemma [2] 
because the determinant Azv (<Zi j <?2 > ' ■ ■ ,Qn-i) given by[3T]has zeros for 8 G (— n, 0). 

Let us consider a case study of a two-site soliton with m = and n<i = m > 2. For clarity of 
presentation, we only consider p > 2. The power series expansions ([7]) give 



m > 3 : 



(<5«,o + Km) (1 + 2e - 2e 2 ) + e 3 W n , n G Z, 



and 



m 

12/ 



2 : 



l2 P 



(52) 
(53) 



{S n ,o + Km) (1 + 2e - 3e 2 ) + e 3 W n , n G Z, 

where {W n } n( zz G / 2 (Z) is a new potential such that ||W||ja = O(l) as e — >• 0. 

Let us consider the coefficient matrix A(9, e) at the continuous spectrum [1, 1 + 4e] defined by 
(f29j) . We have explicitly 



M(5) 



1 e 

-im0 i 



irnO 



N{k) 



1 

-2k 



-2k 



1 



e 

Note that detM(0) = 1 - e _2im9 . Besides the end points 6 = -vr and 6 = 0, the matrix M(6) 
(and, therefore, the limiting matrix A(9, 0)) is singular at the intermediate points #j = — ^ for 
j = 1,2,..., m - 1. 

If m = 2, there is only one intermediate-point singularity of A(9, 0) at 
dimNulU(-|,0) = 1 and 



We have 



NullA* 



7T 



-,0) =span{ei} 



ei 



The first two entries of the right-hand-side vector h(8, e) in the linear system (|30p are given explicitly 
by 

/H(0, e) = J^ e- ie ^fn, h 2 (6, e) = J2 e^^U 



n&L 
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The constraint (ei,h(— ^,0))c4 = of Lemma [3] gives hi{— §,0) = — /i2(— §,0) and it is equivalent 
to the constraint f\ = 0. If / € Z (Z) with f\ ^ 0, then the solution of the linear system (J23H 
and hence the resolvent operator (f22~|) has a singularity at Q, = 1 + 2e (0 = —5) as e — >• 0. This 
singularity indicates a resonance at the mid-point of the continuous spectrum in the anti-continuum 
limit. 

We would like to show that the resonance does not actually occur at the continuous spectrum 
if e > and does not lead to (unstable) eigenvalues off the continuous spectrum. To do so, we use 
the perturbation theory up to the quadratic order in e. 

Expanding solutions of the transcendental equation 



we obtain 



and 



2e(cosh(K) - 1) = 2 + ecv, w = 2 - 2cos(0), 
\ e _^ e 2 + 0(e 3 ) as e ^ 



2 



2esinh(K) = 2 + (2 + w)e - e 2 + C(e 3 ) as e -»■ 0. 
Using expansion (|53p for m = 2, we obtain the extended coefficient matrix A(9, e) in the form 



A(9, e) := 



2ie sin(0) J - (1 + p)v(e)M{9) -pv(e)M(9) 

-pi/(e)JV(«) 2e sinh(«)J - (1 + p)i/(e)JV(«) 



where i/(e) = l + 2e-3e 2 + C(e 3 ). Using MATHEMATICA, we expand roots of dei A(0,e) = near 
9 = — ? and e = to obtain 



7T ... ■> - \2 2 



+ (p-l)e + 2(l-p)e 2 + i(p-l)^e^ + 0(e d ) as e -»• 0. (54) 

Since Im(#) > for small e > and z(X + ) = 6, the solution of the linear system ([30|) is singular at 
the point z(A + ), which does not belong to the domain Imz(A+) < and hence violates the condition 

CZD- 

The singularity of the solution of the linear system fj30|) is still located near the continuous 
spectrum for small e > and, therefore, the resolvent operator R(£l) becomes large near the points 
SI = ±(1 + 2e) (although, it is always a bounded operator from / 2 (Z) x l^(Z) to l^_ a (Z) x ^(Z) for 
small e > and fixed a > 5). Since sin(0) is nonzero for = — ^, the norm of i?(fi) is proportional 
to the 2-norm of inverse matrix yl _1 (0,e). 

Figure [2] illustrates the singularities of the resolvent operator R(£l) by plotting pseudospectra of 
the coefficient matrix A(Q, e) in the complex f2-plane for p = 2 and e = 0.05. The subplots (a) and 
(b) for m = 1 show that the matrix is singular at the edges of the continuous spectrum Vt = ±1 and 
SI = ±(l + 2e), and at four points on the imaginary axis, the latter being attributed to the splitting 
of zero eigenvalue in the anti-continuum limit. The subplots (c) and (d) for m = 2 and m = 3 
respectively show that in addition to singularities at the edges of continuous spectrum there are 
also m — 1 local maxima at its intermediate points. This local maxima correspond to the minima 
of detA(Q,e). We also notice the wedges on the level sets as they cross the continuous spectrum 
occuring due to the jump discontinuities in z(A+) and A(fl, e) because the resolvent operator R(Q) 
is discontinuous across the continuous spectrum. 

Figure [3] further illustrates what exactly happens at the continuous spectrum. On the left, we 
plot L4.(fi, e) L versus 9 G (— 7T, 0) for the case m = 2. On the right, we show that the height of 
the local maxima near 9 = —ir/2 is proportional to e~ 2 as prescribed by formula (JSi 
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Figure U] gives an illustration for pseudospectra of the resolvent operator R(Q). Recall that 
on the continuous spectrum 0, £ [1,1+ 4e], R(fl) is a bounded operator from ^(Z) x / 2 (Z) to 
I'L CT (Z) X /^_ CT (Z)) for fixed <r > 5 ■ To incorporate the weighted I 2 spaces, we consider the renormalized 
resolvent operator 

R L {n) = {L- flh)' 1 : Z 2 (Z) x / 2 (Z) -> Z 2 (Z) x / 2 (Z), 

where L is derived from L by replacing operators /, A and V with /, A and V, and I<i = diag{/, /}. 
Here 

j&u+uu- 

and K n = (1 + w?) u ' 2 . The lattice problem is considered for 2i^ + 1 grid points and the correspond- 
ing matrix representation of operators L and I2 is constructed subject to the Dirichlet boundary 
conditions. 

The level sets for the (2K + 1) x (2K + 1) matrix approximation of the resolvent R(Q) are 
plotted on Figure HI The subplots of Figure [4] correspond to the subplots of Figure [2j We observe 
that the norm of R(Q) has the same global behaviour as for the norm of A(f2,e) -1 . However, the 
resolvent operator R(£l) has no singularities at the edges fi = ±1 and SI = ±(1 + 4e) because these 
singularities are canceled according to Lemma |4] (which remains true for any m > 1, see Remark [6]). 

Although no arguments exist to exclude resonances at the mid-point of the continuous spectrum 
for the linearized dNLS equation ([8|) , the case study of a two-site discrete soliton suggests that the 
resonances do not happen at the continuous spectrum for small but finite values of e > 0. Moreover, 
the resonances do not bifurcate to the isolated eigenvalues off the continuous spectrum because 
isolated eigenvalues near the continuous spectrum would violate the count of unstable eigenvalues 
(I12j) . Therefore, the only scenario for these resonances is to move to the resonant poles on the wrong 
sheets Im(,z(A±)) > of the definition of z(X±). 



A Perturbative arguments for the cubic dNLS equation 

We recall the coefficient matrices A± (e) from the proof of Lemma [21 In the case p = 1 (the cubic 
dNLS equation), these matrices are rewritten in the form 



-2M± -M± 

N(k±) 2esmh(K±)I-2N(K±) 



A ± (e) = 

where k± > are uniquely defined by 

2e(cosh(ft + ) — 1) = 2, 2e(cosh(«; 



l) = 2 + 4e. 



We recall that Null(.A±(e)) and Null(M-i-) are (N — l)-dimensional for any e G [0, eo). It is clear 
from the explicit form of A± (e) that 



weNun(4(e)) ^ u 



w 




w G Null(Mi 



At e = 0, we also recall that Null(A-|-(0)) 2 is (2N — 2)-dimensional because of (JV— 1) eigenvectors 
and (TV — 1) generalized eigenvectors, 



At(0) 




w 



, Afc(0) 



-w 




w 



w 6 Null(M±). 
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m = 1,p = 2, e = 0.05 



m = 1, p = 2, e = 0.05 




0.15 



0.15 



m = 2, p = 2, e = 0.05 





0.15 



E 0- 



(b) 

m = 3, p = 2, E = 0.05 




(c) (d) 

Figure 2: Level sets for L4.(J2, e) L in the fi-plane. The levels are equidistant on a logarithmic 
scale. 



p = 2, £ = 0.05 



300 
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200 
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50 




h=c£ : c = 1 .0671 , a = -2.0048 



-3 -2.5 -2 -1.5 -1 -0.5 



10 




Figure 3: Left: Norm L4(0, e) 1 L versus G (— n, 0) for m = 2. Right: The value of local maxima 
of L4(n, e) _1 || 2 i n the neighborhood of 6 = — rr/2 as a function of e. 
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m = 1,p = 2, <*= 1,6=0.05, K=100 



m = 1,p = 2,(j=1,e = 0.05, K = 1 00 



E 




0.15 



E 




0.15 



E 



m = 2, p = 2, o=1,e = 0.05, K = 100 




0.15 



0.05 



E 



-0.05 



-0.15 



(b) 



m = 3, p = 2, o = 1,6 = 0.05, K= 100 




(c) (d) 

Figure 4: The level sets of (L — O/2 ) — 1 i n the Q- plane. The black dots represent eigenvalues of 
the matrix representation of operator L. The levels are equidistant on a logarithmic scale. 
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We would like to show that Null(At(e)) 2 = Null(At(e)) is (iV-l)-dimensional for any e G (0, eo). 
In other words, we would like to show that no solution u G C of the inhomogeneous equation 
A±(e)u = u G Null(^4-|-(e)) exists for e G (0, eo)- This task is achieved by the perturbation theory. 
We will only consider the case A + (e), which corresponds to 6 = 0. The case A__(e) which corresponds 
to 6 = —it can be considered similarly. 

We shall only consider the case of the simply-connected set C/+ U U- with m\ = mi = ... = 
771 A?- 1 = 1- The general case holds without any changes. 

Thanks to the asymptotic expansions 



-K+ 



+ 0(e 2 ), 2esinh(K + ) = 2 + 2e + C(e 2 ), 



as 



0, 



we obtain the asymptotic expansion 
A + (e) = 



2M 4 
-I 



-M 4 
O 



+ e 



O 



O 



I J 21 -J 



+ OU 2 



where I and O are identity and zero matrices in R and J is the three-diagonal matrix in 



J 



>N 






1 


• 


• 





1 





1 • 


• 








1 


• 


• 











• 


• 


1 








• 


• 1 






(55) 



Note that (2/ — J) is a strictly positive matrix because it appears in the finite-difference approxi- 
mation of the differential operator — d% subject to the Dirichlet boundary conditions. 

Perturbative computations show that if u G Null(^4+(e)), then u is represented asymptotically 

as 



e(2/ - J)v 



+ OU< 



where v + 2e(2J - J)v + C(e 2 ) = w G Null(M+). 

Now, there exists a solution u G C 2N of the inhomogeneous equation A + {e)u = u£ Null(A + (e)) 
if and only if u _L Null(^4^_(e)). For small e G (0, eo), this condition implies that 

e(2/ - J)v + C(e 2 ) = e(2/ - J)w + C(e 2 ) J_ w G Null(M+), 

which is not possible since (21 — J) is a strictly positive matrix. 
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